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; ^ , Abstract. It is shown that the ground state energy of heavy atoms is, to 

■ leading order, given by the non-relativistic Thomas-Fermi energy. The proof is 

based on the relativistic Hamiltonian of Brown and Ravenhall which is derived 
from quantum electrodynamics yielding energy levels correctly up to order 

(N 

> 

Q\ ' 1. Introduction 

l> 

, The energy of heavy atoms has attracted considerable interest in the context of 

CO ' nonrelativistic quantum mechanics. Lieb and Simon |2(J| proved that the leading 

behavior of the ground state energy is given by the Thomas-Fermi energy which 
decreases as Z'^^^. The leading correction to this behavior, the so called Scott 
, correction was established by Hughes El (lower bound), and Siedentop and 

' Weikard O EHl EEl 123 EHI (lower and upper bound). In fact even the existence 

of the ^5/3 -correction conjectured by Schwinger was proven (Fefferman and Seco 
' [51 IIOL rm 1^ 1121 171 151 1^ 15]). Later these results where extended in various ways, 

, e.g., to ions and molecules. 

Nevertheless, from a physical point of view, these considerations are questionable, 
since large atoms force the innermost electrons on orbits that are close to the nucleus 
where the electrons move with high speed which requires a relativistic treatment. 
, Our main goal in this paper is to show that the leading energy contribution is 

' unaffected by relativistic effects, i.e., the asymptotic results of Lieb and Simon pTI) 

remain also valid in the relativistic context, whereas the question mark behind the 
quantitative correctness of the other corrections persists. 

S0rensen |23j took a first step in this direction. He considered the Chandrasekhar 
multi-particle operator and showed that the leading energy behavior is given by the 
non-relativistic Thomas-Fermi energy in the limit of large Z and large velocity of 
light c. Nevertheless, a question from the physical point of view remains: Although 
the Chandrasekhar model is believed to represent some qualitative features of rela- 
tivistic systems, there is no reason to assume that it should give quantitative correct 
results. Therefore, to obtain not only qualitatively correct results it is interesting, 
in fact mandatory, to consider a Hamiltonian which - as the one by Brown and 
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Ravenhall - is derived from QED such that it yields the leading relativistic 
effects in a quantitave correct manner. 

2. Definition of the Model 

Brown and Ravenhall |2] describe two relativistic electrons interacting with an 
external potential. The model has an obvious generalization to the A^-electron case. 
The energy in the state -ip is defined as 

N 

(1) "^^^ 

u=l l<fi<u<N 

where 

Dc,z :=«• yV + c2/3-Z|-|-i 

is the Dirac operator of an electron in the field of a nucleus of charge Z. As 
usual, the four matrices ai, ...,Q!3 and /3 are the four Dirac matrices in standard 
representation. We are interested in the restriction £ of this functional onto Qn ■— 
A^^i(ff^/^(K^) ® C'') n S^N where 

N 

(2) SjN ■■= A ^' 

the underlying one-particle Hilbert space is 

(3) ii:=[x(o,oo)Pc,o)](i'(R')®C4). 

Note that we are using atomic units in this paper, i.e., me — h — e — \. 

As an immediate consequence of the work of Evans et al. '3' this form is bounded 
from below, in fact it is positive (Tix |29ll3(Jj ). if k Z/c < K^iit '■= 2/(7r/2 + 2/7r). 
(In the following, we will assume that the ratio n e [0, Kcrit) is fixed.) According 
to Friedrichs this allows us to define a self-adjoint operator Bc,n,z whose ground 
state energy 

(4) E{c, N, Z) := inia{B,,N,z) = inf{£(V')|V e Qn, U\\ = 1} 

is of concern to us in this paper. In fact - denoting by Etp{Z, Z) the Thomas-Fermi 
energy of Z electrons in the field of nucleus with atomic number Z and q = 2 spin 
states per electron (see Equations (flTj) and (fTH|l for more details) - our main result 
is 

Theorem 1. 

E{Z/k, Z, Z) = Etf{Z, Z) + o{Z'^'^). 

This result, given here for the neutral atomic case, has obvious generalizations 
to ions and molecules. To keep the presentation short we refrain from presenting 
them here, as their treatment follows the same strategy. 

The remaining paper is structured as follows: First we show how the treatment 
of the Brown-Ravenhall model can be reduced from Dirac spinor (4-spinors) to 
Pauli spinors (2-spinors). In Section |3| we prove the upper bound corresponding 
to Theorem n] by rolling it back to Lieb's upper bound in the non-relativistic case 
[T7|. Section ^reduces the lower bound to S0rensen's lower bound Finally, in 
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the appendix we show that the correlation estimate using the exchange hole yields 
a pointwise lower bound with uniform error of order Z . This is interesting in itself 
since it allows to estimate the error purely by the particle number not using any 
kinetic energy. 

We now indicate, how to reduce to Pauli spinors. To this end we parameterize 
the allowed states: Any ip G Sj can be written as 

(5) W, ] 

for some u E t) := L^(R"^) (g) C^. Here, cr are the three Pauli matrices, 

p := -iV, E,{p) (c2p2 + c^^^, iV,(p) := [2E,{p){E,{p) + c^)]"\ 
In fact, the map 

$ : [} ^ i3 




(6) (-^lu 

embeds t) unitarily into and its restriction onto H^{E.^) ® €? is also unitary 
mapping to^)C^ ^^(R^) ^ (Evans et al. HJ) ). 
It suffices to study the energy as function of u 

N 

(7) £ o («li<i>) : /\ b -> R. 

i/=i 

The one-particle Brown- Ravenhall operator By for an electron the external elec- 
tric potential of a point nucleus acting on Pauli spinors is then 

(8) B,^z-^E,{p)-Z^i-Z^2 
where we have split the potential into 

(9) (pi :=$J|-ri$i, (^2:=*;i-r'$2. 

As we will see the first part ifi is contributing to the nonrelativistic limit whereas 
the second part turns out to give energy contribution that do not even affect the 
first correction term. 

3. Upper Bound 

3.1. Coherent States. The upper bound will be given by choosing a trial density 
matrix in the Hartree-Fock functional for the Brown-Ravenhall operator. To this 
end we introduce spinor valued coherent states. 

Given any function / e iJ^/^(R^) and an element a — (p, q, r) of the phase 
space r := R^ X R'^ X {1, 2}, we define coherent states in t) as 

(10) Fa{x) := (v?p,q «) e^)(a;) := /(x - q) exp(jp • x)(5r,^, 

where x — (x, cr) e R"^ x {1,2} and the vectors Ct are the canonical basis vectors 
in (see Lieb ^7] and Evans et al. 3 ). We will pick / depending on a dilation 
parameter. More specifically we will choose 

(11) /(x) to(x) := i?-3/2g(i?-ix) 

where R :— Z^^ with 5 e (1/3,2/3) and g € H^^^, spherically symmetric, normal- 
ized, and with support in the unit ball. 
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The natural measure on F counting the number of electrons per phase space 



volmnc in the spirit of Planck is Jp dn{a) :— {2n) / dp / dq J2t=i- "^^^ essential 
properties needed are the following. For A e L^{T,dn) 

(12) -f:= J^dn{a)A{a)\F^){F^\ 
is a trace class operator and 

(13) 0<^<1 =^ 0<7<1 

(14) tr7 = J dn{a)A{a). 

Using $ we can lift any such operator 7 to an operator on Sj 

(15) 7$ := $7$*. 
We will pick 

(16) A{a) X{($,x)eK6|^2/2-Vz(x)<o}(Pi q) 

where Vz Z / \ ■ \ — \ ■ \ * Ptf'i here pxF is the unique minimizer of the Thomas- 
Fermi functional 



(17) £tf{p) := / 

Jr3 



^7TFp(x)5/3 _ Z ,^ 



dx + D{p, p) 



where, for Fermions with q spin states per particle, 7tf '■— {6n'^ /q^^'^h'^ /{2m), i.e., 
in our units, 7tf = (377^)^/^/2. Note that J dQ{a)A{a) = Z (Lieb and Simon 
mi). Note also that T4(q) := Z'^/^ViiZ'^^^q) (see also Gombas (1^ and [HI). Note 
also that the minimal energy Etf {N, Z) fulfills the scaling relation 

(18) ETFiN, Z) = Etf{N/Z, 1)Z^/^ 

Note, that we could restrict the minimization to J p < N without any problem. 
For N > Z there would be no change in the minimizer; for N < Z we would get a 
different minimizer. For notational convenience we will merely consider the neutral 
case N = Z in the following. 

3.2. Upper Bound. We begin by noting that the Hartree-Fock functional - with 
or without exchange energy - bounds E{c, iV, Z) from above. To be exact we 
introduce the set of density matrices 

(19) Sn := {7 e &\^) \Ec{ph e ©^(f)), < 7 < 1, tr7 = iV} 
where denotes the trace class operators on f). 

ffjF ■ Sn — !■ K 

^^^^ 7 ^ tr[(£;,(p) - c2 - Z/|x|)7<,] + D{p^,,p^,) 



where - as usual - p-y is the density associated to 7 and D is the Coulomb scalar 
product. By the analogon of Lieb's result |18l [T?)j (see also Bach [Jj) - which 
trivially transcribes from the Schrodinger setting to the present one - we have for 
all J e Sn 

(21) E{c,N,Z) <£^F{-f). 
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3.2.1. Kinetic Energy. By concavity we have 

(22) E,{i;y) - < ip^ 

which imphes that the Brown-RavenhaU kinetic energy is bounded by the non- 
relativistic one, i.e., for all 7 G S'jv with — A7 e S^([)) 

(23) tr[(i?,(p)-c2)7] <tr(-iA7). 

Inserting our choice of 7 (see Equations lfTU|l . ((TT|l . ((T^ . and turns the 

right hand side into the Thomas-Fermi kinetic energy modulo the positive error 
Z||V5||2i?-2 (ggg LjgT^ [Y7| Formula (5.9)]), i.e., 

(24) tr[(£;,(p) - c2)7] < ^7tf j Ptf Wdx + ZR-^\\Vg\\\ 

3.2.2. External Potential. Since — Ztr((/327) is negative, we can and will estimate 
this term by zero. This estimate will be good, if this term is of smaller order. 
Although, logically unnecessary for the upper bound, it is, for pedagogical reasons, 
interesting to see that if2 does indeed not significantly contribute to the energy, if 7 
is chosen as above. Moreover, the proof will be also useful for the proof of Lemma 

m 

Lemma 1. For our choice of j — J-^dfl{a)\Ea){Fa\ and 6 e (1/3,2/3) we have 
(25) 



< Ztr((^27) <kZ dn{a)A{a) 



(In the following - throughout the paper - we use the letter k for a constant 
independent of c, N, R, or Z.) 

Proof. We begin by estimating the expectation of ip2 in a coherent state. 



(26) 0< (F„,(p2i^a) <A: // d|d| 



< kc-^R-' 1 1 d^d^'^^^^^\^ + RpW + Rp\ < kc-^R-Hl + R\p\ + R'\p\^). 



..1.9(0.9(^)1 

Here, we used that Nc{^) > and, in the last step, that 

is integrable in ^ and ^' because g G i/'^/^(R'^). Thus we get 



(27) < Ztr((^27) - d<^{a)A{a)(F^,^2Fo) 

<fc^ / d"f^(«M(«)(i+^ipi+^'ipn 

<A;-^|z + i? I dq [z4/3v^i(Zi/3q)]' + i?2 \^Z^I^Y^{Z^l^c^'^^ 



□ 
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Lemma 2. For our choice of "f and 6 € (1/3,2/3) we have 
(28) |Ztr[(|.|-^-^i)7]| 

<-i.7 [AO(r.^A(r.\ [[ ^^^(^ (^c(0+C^)(igc(O+C^) VA ...nA 

= 0(^5/3+*). 

Proof. We fisrt note that 



(29) 



^ (i^e(g)+c^)(£,(€0+c^) 

iVc(^)iVe(r) 



< 



3ii;e(4)i?c(r) - c'iEciS,) + £^0(4)) + C2 



Nc{i)Nc{i') 

Then, noting that Ec{£,) ~ < c|^|. we obtain 



(30) 



1 - 



Nc{i)Nc{i') 



< 



3c^|^|iri+2c3(^ + 



2c4 

Using this last equation, we estimate 
(31) |(F„,(^-^i)J^„)| 

-'JJ V mm^) — ) "^"^^^""^"(^ ^" 



<fc / dCd^'Mi^^)^|^(e-|^||^'|+e-(|^| + |C'|)) 

<kc-'R-'j y d^'k^M^d^+ijpii^'+ijpi+ci^i^+ijpi) 



<kc-'R-'J 1 d^'^^M^(|||||'|+i?|p||^| + |i?p|2 + ci?|||+ci?2|p|) 

< fcc-2i?-3(l + R\p\ + i?2|p|2 + ci? + Ci?2|p|). 

Thus 

(32) Z\ tr[(| • |-i - ^1)7] < Z| y d-Jl(a)A(a)(F„, (I • \-' - ^i)F„)| 

< fc(Z3^ + ^2^+2/3 + +^5+4/3 + ^25 ^ ^5+5/3) 

which yields the desired estimate. □ 

3.2.3. The Electron- Electron Interaction. We will roll back the treatment of the 
electron-electron interaction to the treatment of nucleus-electron interaction. 

Lemma 3. For our choice ofj and 5 G (1/3,2/3) we have 

(33) £'(P7*. P7*) - D{p^, Pi) = 0{Z'>^^+'), 



RELATIVISTIC ENERGY 



7 



where pj is the density of 7 and p-y^ is the density of 7$ . 
Proof. We have 

(34) \mP'r+pi.)*\-r'm\ < Vv^wpy+p^jm-' 

Now, 

\D{p.y^,p-f^) - D{p-^,p^)\ ^ \D{p-^^ ~ Pl,P-l^ + P'l)\ 
1 

< - 



23/2^-l/2^|^|- 



(P7 W - W)[(P7 + P7*) * I • I ^](x)dx 

<\ j^dn{a)A{a) 
1 1 d|dr|^[(p, +p^J * I • |-i](| -|')l^(ei')|i^a(l)l|i^a(l')|d|d|' 
< j^dn{a)A{a) j j didi'U- i'\-^K{i,£:)\F^{i)\\F^{i')\didi' 



where 



{E,{i) + c^){E,{i')+^) 



and where we used (|34|) in the last step. Eventually, applying Lemmata ^ and |21 
yields the desired result. □ 

3.2.4. The Total Energy. Gathering our above estimates allows us to reduce the 
problem to the non-relativistic result of Lieb T7! 

Theorem 2. We have E{kZ, Z, Z) < ETpil, 1)^^/^ + kZ^°/^. 

Proof. Following Lieb jl7l Section V.A.I] with the remainder terms given there 
(putting R = Z^^ as in our estimate), using the remainder terms obtained in 
Lemmata n through 131 and using l|24|l we get 

(35) E{c, Z, Z) < £hf(7) < Etf{Z, Z) + 0{Z^+^^ + Z^'i + zi+^) 
which is optimized for (5 = 5/9 giving the claimed result. 

4. Lower Bound 



□ 



The lower bound is ~ contrary to the usual folklore - easy. As we will see, it is a 
corollary of S0rensen's [221 result for the Chandrasekhar operator and an estimate 
on the potential generated by the exchange hole [2Tj . The exchange hole of a density 
cr at a point x g K'^ is defined as the ball i3fl;^(x)(x) of radius Ra{^) centered at x 
where Ra{^) is the smallest radius R fulfilling 

(36) i = 



The hole potential of a is defined through 



(37) 



i.(x) := / 



^(y) 



/|x-y|<J?,„(x) |X-y| 

Our second main result is the following lower bound. 



■dy. 



8 



R. CASSANAS AND H. SIEDENTOP 



Theorem 3. 

liminf[(£;(c, Z, Z) - Etf{Z, Z)\Z-''/'^ > 0. 

Z^oc 

Proof. Pick 5 > and set ps ptf * 9^-^- Then the exchange hole correlation 
bound (211 Equation (14)] implies the following pointwise estimate 



1 

Because of the spherical symmetry of g we can use Newton's theorem and 
replace ps by ptf in the third summand of the right hand side of H38I) . Then, by 
Lemma |S1 we get that for all normalized Cz Qn 

(39) f (^) > tr[A+(|i?o| - - Vs)A+]- - kNZ - i?(pTF, Ptf) 

where, for t G M, we set := min{t, 0} and Vg — Z/\ ■ \ — ps * \ ■ \^^. 

To count the number of spin states per electron correctly, i.e., two instead of the 
apparent four, we use an observation by Lieb et al. Appendix B]: Note that 



(40) A_ = U-^A+ U, where U := 

Indeed, we have 



1 
-1 



and 



We set X := {\Do\ - - Vs{x))l2, and write 
tr 








a:_ 


))- 




- 


(o 


a:_ 


= tr ( 


A.( 



Thus 

(41) 2tr(A+ 



Since ji^ol = Ec{p), we obtain 

(42) EitiZ, Z,Z) >2 tr[£;,(p) - - V^5(x)]_ - i^(pTF, Ptf) - kNZ 

where the last trace is spinless. This connects to S0rensen's Equation (3.2) from 
[221 ■ The result then follows using his lower bound. □ 
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Appendix A. L°°-Bound on the Exchange Hole Potential 

We begin the appendix with the following remark: the Thomas-Fermi potential 
Vz ■= Z/\ ■ I — ptp * I • |~"^ can be written as 

(43) 7tf4f = yz 

(see, e.g., Gombas This equation yields immediately the upper bound 

(44) ptf(x) < (Z/7tf)3/2|x|-3/2. 

This bound allows us to prove the following L°°-bounds on potentials of exchange 
holes. 

Lemma 4. 

||L,,JU = 0(Z). 

Proof. The function 

/ : M+ ^ K 

(45) 



Vt f |yri|y+(0,0,l)|-3/2dy 

J\y\<l/t 



'|y|<i/t 

is obviously continuous on (0,cx)). Moreover, f(t) tends to a positive constant for 
t ^ and to for t oo. Thus, ||/||oo < oo. 
This allows us to obtain the desired estimate: 

(46) Lp^p(x) < Ai(x) + A2(x) 

where 

(4^) ^^w-/ ^^dy<m''7 I II 1 13/. 

J\y\<i/z \y\ V7TF/ 7|j,i<i/z |y||y + x|3/2 

= {Z/jTF)'/'Z''/'f{\^\Z) < ||/||oo7^F^'^. 

and 
(48) 

PTF(x + y) , ^ „ f r . \j ^ ^ 



A2{y) ■■= / I M y < Z / ptf(x + y)dy < . 

"'^<|y|<-R.PTFW 1^1 "'^<|y|<-RpTFW ^ 

These two estimate proof the claim. □ 

Lemma 0] allow us already to estimate the N electron operator Bc,n,z by the 
canonical one particle Brown-Ravenhall operator whose nuclear charge is screened 
by the the Thomas-Fermi potential. However, since we would like - because of mere 
convenience - to take advantage of S0rensen's result we derive an estimate on 
Lp_5 (where ps '■— ptf * dz-s), i-e., the exchange hole potential of the density 
occurring in S0rensen's proof. 

Lemma 5. 

||ip,IU = o(z). 
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Proof. We proceed analogously to the proof of Lemma 01 

P5(x + y), , f ps{x + y) 



(49) L,.(x)< / ^^i^ftZldy+Z 

J\y\<l/Z \y\ Jl 



1/Z<|y|<i?p (x) |y| 



dy 



< / dz5|-.(z) / Pll(^l_l±ylay + Z f psix + y)dy 

< J dz.g|_, {z)Ai {x- z) + ^ < kZ 

where we used the definition of the radius of the exchange hole from first line to 
the second line, the definition of Ai in the next step, and in the last step the 
L°°-estimate on ^i. □ 
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